Producing GHZ state of nitrogen-vacancy centers in cavity QED 
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We propose a scheme to generate GHZ state of nitrogen-vacancy centers coupled to a whispering-gallery 
mode cavity. The implemented evolution is independent on the cavity field state with the assistant of a strong 
classical field, and thus not sensitive to the thermal state of the cavity. Meanwhile, it is fast compare to the 
convectional dispersive interaction induced operation in a cavity-assisted system. The scheme is readily scalable 
to multiqubit scenario. 
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Entanglement is one of the most counterintuitive proper- 
ties of quantum mechanics, which plays an central role in 
quantum information and quantum computation 1 1]. For real 
quantum system, it is generally sensitive to practical noise, 
which results from the infamous decoherence effect. There- 
fore, generation of entangled state with inherent robustness to 
diverse noise is preferred. Due to long electronic spin lifetime, 
fast initialization, good qubit readout, and coherent manipula- 
tion at room temperature, the diamond nitrogen-vacancy (NV) 
center is considered as a promising candidates for quantum 
computing |2]. Coupling of many qubits in solid state sys- 
tems is problematic. For NV center, the partial solution is in- 
troduce coherent coupling between electron and nuclear spin 
qubits. However, this can only be done in the case of a few 
qubits since there are limited number of nuclear spins can be 
individually addressed [3]. Alternatively, separated NV cen- 
ters can be coupled by a cavity mode |4]. For a easy fabricated 
microsphere cavity, the quantized whispering-gallery mode 
(WGM) may has a ultrahigh Q factor and very small volume 
list], and thus leading to reach the strong coupling regime. En- 
tanglement generation schemes are also proposed with a mi- 
crosphere cavity |6] or a one dimensional resonator [7]. How- 
ever, the previous schemes all work in the dispersive regime, 
i.e., the detuning between the cavity mode frequency and the 
qubits is much larger than that of the cavity-qubit interaction, 
which makes the effective coupling relatively weaker com- 
pare to resonant coupling, that is, the generation time will 
be longer Meanwhile, the coupling usually contains cavity- 
state-dependent energy shift, which leads the generation to be 
sensitive to the cavity thermal state. 

Here, we propose an alternative scheme for generating mul- 
tipartite entangled state of nitrogen-vacancy centers coupled 
to a whispering-gallery mode cavity, which overcome the 
above-mentioned drawbacks. In addition to the optical transi- 
tions driven by the cavity mode and optical laser, we introduce 
another transition within the ground states strongly driven by 
a microwave field. Before proceeding, we want to emphasize 
that using both optical and microwave radiation to control an 
electron spin associated with the NV center in diamond is ex- 
perimentally accessible |J8J] . With the driven of a microwave 
field, one get two additional merits of the proposed scheme. 



Firstly, the photon-number-dependent parts in the evolution 
operator are canceled, and thus the scheme is insensitive to the 
thermal state of the cavity mode. Secondly, the large detuning 
constrain can be removed so that fast gate operation can be 
achieved. Meanwhile, the solid-state set-up is readily scalable 
to multiqubit scenario and the gate speed is not slowed down 
with the increasing of the involved qubits. 

As illustrated in Fig. 1(a), A^ identical NV centers in N 
separate diamond nanocrystals may be strongly coupled to 
the WGM of a microsphere cavity. The lowest-order WGM 
corresponding to the light traveling around the equator of the 
microsphere, which offers exceptional mode properties for 
reaching strong coupling regime. One can model the NV cen- 
ter as a three-level system, as shown in Fig. 1(b), where the 
states pA, TTig = O) and ^AjTO^ = — l) are encoded as our 
qubit states |0) and |1), respectively. The state \^E,ms = O) 
is labeled by |e) and do not use the metastable ^A state, which 
has not yet been fully understood H. In our implementation, 
the optical transition |0) —>■ |e) and |1) — > |e) (with transi- 
tion frequencies ujeo and ^ei ) are coupled by the WGM with 
frequency cUc and a laser with frequency ujl and polarization 
(T+ 1 10], respectively. Both coupling is far-off resonant from 
their transition frequencies so that the |e) state can be adia- 
batically eliminated. The NV centers are fixed and apart with 
the distance much larger than the wavelength of the WGM, 
so that they can interact individually with the laser beams and 
the direct coupling among them is negligible. Then, in units 
of h = 1, the system is desciibed by 



Hs = ^ca+a + y^a;»|z)(i| 



(1) 



iV 



+ ^[G,a|e)(0| + f7L,,e-^("-+'^)|e)(l|+H.c.], 

where a'^{a) is the creation (annihilation) operator of the 
WGM field, oji is the frequency of ith energy level with 
i € {0, 1, e}, Gj and flL,j is the coupling strength between 
jth NV center and WGM and laser, respectively. 

Using the rotating-wave approximation (RWA), with re- 
spect to the free Hamiltonian Hq = ujcO^a + ^iOJi\i){i\, 
the interaction Hamiltonian in the interaction picture can be 
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m =0 
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where a^ denotes Pauli matrix in the eigen basis of (Jx with 
eigen equations as ax|±) — ±|±), where |±) = (|0) ± 
|l))/\/2- In the case of 17 ^ {<5, ??}, we can omitting the 
fast oscillation terms with frequencies Q,±6 using RWA, then 
the Hamiltonian (JUi reduces to 

i73=r;(ae-^(^'*-^)+aV(^*-^))j., (5) 

where J^ ~ '}2j=i "'? /2- For = 0, the above Hamiltonian 
reduces to 



FIG. 1 : (a) Schematic setup of the microsphere cavity, where A'^ iden- 
tical NV centers (red dots) in diamond nanocrystals are equi distantly 
attached around the equator of a single fused-silica microsphere cav- 
ity, (b) Level diagram for a NV center, with G and f^Lare the cou- 
pling strength between NV center and WGM and the laser pulse, re- 
spectively. We encode qubits in the subspace spanned by the state of 
nis = and m^ = —1. The qubit zero-field splitting is tJio =2.88 
GHz and can be driven by a microwave field with strength Q,. 



written as Hill] 
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7J, = ^7?,(aa+e-''(^^*-«+H.c.) 



(2) 



j=i 



where a+ = |1,) (0,|, nj = |0,) (1 
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-Aj+Sj ' . 

ijjc — i^io,j — '^L.j- For simplicity, hereafter, we assume that 
r\j — rj and Sj — S. However, this is not necessary here. The 
inhomogeneous qubit-cavity coupling will result in deviation 
of the generated final state from the target state by extra phase 
factor, which can be easily compensated by single qubit op- 
erations lll2ll . Meanwhile, besides the optical coupling, we 
also add a resonant microwave driven field at a frequency of 
the microwave transition ojiq for each NV center Then, in 
the above interaction picture, it reduce to Hd — ^c^ with 
the Pauli matrix written in the qubit subspace. Then the total 
interaction Hamiltonian can be written as 



H, = 



N 

E 






aa+e-'(**-^) + 



a+a^ e 



i{St-(t>) 



(3) 

This microwave control over the NV centers is well studied 
experimentally [13] and even all the three magnetic sub-levels 
can be resolved 1.141 . With this strong microwave coupling, 
we can realize the analogy of str ong -driven atomic cavity 
QED system with two level atoms [15]. IVIean while, this dis- 
tinct our scheme with all the previous investigation with NV 
centers coupled by a cavity mode [IJ, la lZ|| . In the interaction 
picture with respect to the first term, the interacting Hamilto- 
nian reads lll5l] 
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The time-evolution operator for the Hamiltonian in Eq. (|6} 
can be expressed as lll6l - ll9[] 

U{t) = exp [-i^(i) J^] 

X exp [-'iB{t)aJx\ exp \—iB*{t)a)Jx\ , (7) 

where 



Ait) 



T 



1 

Ts 



iSt 



-l)-t 



Bit) 






-iSt 



1 



(8) 



(9) 



While the whole time evolution is nonperiodic, Eq. (|9j shows 
that Bit) is a periodic function of time, and it vanishes at 
5T = 2fc7r where k = 1,2,3,.... At these times, the time 
evolution operator in the interaction picture reduces to 



[/(Y) = exp(z7' J2), 



(10) 



with i = TfT/S = XT. 

The operator in Eq. (fTOt is readily to be used for gen- 
erating GHZ state. However, this requires a carefully cho- 
sen duration of the interaction because it involves real ex- 
citation of the cavity photon. Fortunately, idea similar to 
the spin echo and dynamical decoupling can be used to op- 
timize the generation process [ 20. - 22] . The key is the map of 
Hdi -^ H£i2 = —Hdi, which can be conveniently achieved 
by setting = tt in our scheme. Then, a gate 



Uil) = exp(i7J^), 



(11) 



with 7 = 'V ( I ^^^ T ^ ^) ^^^ ^^ obtained as the following 
1 20I - I22I] . For a gate with duration t, one can set Hjji for the 
duration of '-^ t/2 and Hd2 for the left duration t/2 ^ t 
1201] . similar to the dynamic decoupling method. Alterna- 
tively, the map of Hdi -^ —Hdi can be done by it pulses 
on the qubits 1,121 . similar to the spin echo method. 

It is notable that C/(7) operation proposed here possess the 
following distinct merits. Firstly, it is insensitive to the ther- 
mal state of the cavity mode as the related influence repre- 
sented by the last two exponents in Eq. (|7|l is completely re- 
moved. Secondly, it also remove the constrain of large detun- 
ing. It is noted that T ^ 1/ry for S ^ t], which is comparable 



to the resonant coupling strategy and is much faster than that 
of the conventional dispersive coupling with 5 ^ ij. This fast 
gate operation is very important in view of limited coherence 
times. Thirdly, it is obvious that the time needed for this gate 
operation is comparable to a two-qubit gate as A{t) is inde- 
pendent on N, and thus the gate speed is not slowed down 
with the increasing of qubits. Therefore, it is readily scalable 
to multiqubit scenario. 

Now, we illustrate how to generate the GHZ state in this 
system, which has many application in quantum communi- 
cation [.23.1 . By choosing 7 = tt/2 and an initial state 
|*)i = |00 • • • 0) for iV NV centers, the final state 
exp (i-| Jj) |\I')i is found to be a GHZ state given by | 

1^^) = ^[e^*T|00 • • • 0) + e*"(3+f )|ii . . . 1)], (12) 
v2 

when N is even. Note that 7 ~ tt/2 can be achieved when 
choosing the parameter as (5 = 2\/krj. Then, for fc = 1, one 
obtains 5 — 2rj and the entanglement generation time as T = 
■n/rj. For odd N, one can get GHZ state by applying another 
unitary operator U = exp(— i-| J^;) in addition to Eq. ( fTTT i. 

In the derivation of effective Hamiltonian in Eq. (|5]l, we 
have neglected the following terms 
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= yae^*'^* ^ [sin(m)CT| - 008(17^)^}] + H.c, (13) 
i=i 

which could reduce the fidelity of the generated state. How- 
ever, for az terms, one can use the spin-echo technique to 
eliminate the infidelity fl2']. Following Ref. [19], we con- 
sider the effect of the (jy terms on the gate operation. In 
the interaction picture with respect to our free Hamiltonian 
Hfree = Sa'^a + X]i=i T^'^T' ^^^ interaction Hamiltonian is 
Hnj{t) = W{t)HnU{t) with U{t) = exp{~iHfreet), and 
we have the propagator U'{t) from the Dyson series. During 
the integration, one can treat U{t) as a constant as Hn{t) is 
oscillating much faster than the propagator Then, under this 
assumption, we get 
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2n 

2 '^ 



(14) 
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^[1 - cosi2nt)]UHt)a^alU{t)} + ■ 
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At the time 5t — 2mr, the infidelity is 

F,n = ^[l - cos{2nt)], 



(15) 



where £, = N{N - l)ri'^ / {8fl'^) . Noted that the estimation of 
Fin is obtained in the interaction picture. Return to the rotat- 
ing frame, the time evolution operator will get the following 



additional term of C/^; — exp(— i^ ■ j^ ^cr^^i). Therefore, 
control the duration of classical microwave field accurately to 
fulfill iU = AriTT, the effect of the a-y terms and Ux will both 
vanishes. Note that 5t = 2mr, so ft/S — 2n. Considering the 
further requirement ofVt^d, we may choose Q, — 6(5. 

To see the feasibility of the present scheme with current 
technology, we now use typical values of physical parame- 
ters in the system to estimate the operation time T. NV cen- 
ters are located near the microcavity surface in order make 
the coupling between the NV center and WGM maximum, 
which could be Gmax — 27r x 1 GHz with the transition wave- 
length between the states |e) being 637 nm, the spontaneous 
decay rate of the excited state being Fg = 27r x 83 MHz, 
and the mode volume be Vm =100 ^im? ll2411 . For (5 ^ A, 
?7 ~ 2GnL/A = 27r X 50 MHz, with f^L = 27rx 500 MHz 
and A = 27rx 20 GHz. The effective qubit spontaneous emis- 
sion rate is estimated as fl^ T^ff = FoJIlG/A^ w 27r x 0.1 
MHz. In the fused-silica microsphere cavity, the small radius 
of 10 /im could lead to the Q factor exceeding 10^ with pos- 
sible improvement as demonstrated in |26|. The cavity decay 
rate rate is k = uJc/Q = 27r x 0.47 MHz for Q = 10^. In our 
scheme, for & ^ 2rj = 2tt x 100 MHz, the assumption 5 <^Q. 
with f7 = 27r X 600 MHz 1 13] and (5 < A is well-satisfied, 
then the entangled state generation time is T = 10 ns, which 
is negligible small compare to both the cavity decay and qubit 
coherence time even at room temperature. 

With the above parameters, for N — A, ^ k, 0.01, and 
thus Fin < 2.1%, which increases with the increasing of A^. 
To see the time dependence of the fidelity of the generated 
entangled state, as shown in Fig. (2), we numerically simu- 
lated the time-evolution of the effective Hamiltonian of Eq. 




FIG. 2: Fidelity of the entanglement generation, see the main text for 
the chosen parameters. 



(11) plus Hn in Eq. (14). We have defined the fidelity as 
F = {^>^U{'j)\'^^) X (1 - F^n). From the lower one of Fig. 
(2), we can see that the generation is now very robust to small 
deviation of the time duration. 

In summary, we have proposed a scheme to generate GHZ 
state of nitrogen-vacancy centers coupled to a whispering- 
gallery mode cavity. The fast scheme is not sensitive to both 
the spontaneous emission decay of the NV centers and ther- 
mal state of the cavity. Finally, the scheme is readily scalable 
to multiqubit scenario as the generation time is independent 
on the number of involved qubits. 
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